Multifragmentation of non-spherical nuclei 
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Influence of the shape of thermalized source on various characteristics of multifragmentation process 
as well as its interplay with effects of the angular momentum and the collective expansion are 
studied for the first time and the most pertinent variables are proposed. The analysis is based on 
the extension of the statistical microcanonical multifragmentation model. 
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In studying the multifragmentation process, a large 
range of incident energies, changing by about four or- 
ders in magnitude, was covered and various types of 
projectiles, from proton till heaviest available ions, were 
probed. The reaction mechanism is often considered in 
terms of two-step scenario where the first, off-equilibrium 
and dynamical step results in the formation of thermal- 
ized source which then, in the second step, decays sta- 
tistically into light particles and intermediate-mass frag- 
ments (IMF's). Assuming that the thermal equilibrium 
is attained, various statistical multifragmentation mod- 
els were employed for the second step (see |l|-|3| and 
references quoted therein). These models were so suc- 
cessful in providing an understanding of basic aspects of 
the multifragmentation process that the deviations be- 
tween their predictions and the experimental data have 
been often taken as an indication for dynamical effects 
in the multifragmentation. This kind of simplistic 'cause 
- effect' interpretation may however be misleading due 
to several oversimplifying assumptions in the statistical 
calculations, such as , e.g., the spherical shape of the 
thermalized source. Indeed, one expects that the spheri- 
cal shape can be perturbed during the dynamical phase 
and the density evolution has both compressed and rar- 
efied zones which can give rise to a rather complicated 
source forms [§J^]. Perhaps more important are the an- 
gular momentum induced shape instabilities |^, 7) which 
may cause large fluctuations of both the Coulomb barrier 
and the surface energy even for moderately high angular 
momenta (L ~ AQTi). Moreover, at high excitations, not 
only the quadrupole stiffness becomes small but also the 
fission saddle point moves towards larger elongations and 
smaller neck cross-sections J?]], giving rise to some 'neck 
effects' ^,0. Hence, before discussing dynamical effects 
in the multifragmentation decay, one should study the 
effects of different shapes in the freeze-out configuration. 
In this paper, the non-spherical fragmenting source is 
considered within the statistical model and the observ- 
ables sensitive to the shape of this source are searched 
for. 

Our statistical consideration is based on the MMMC 
method of the Berlin group Q . In the MMMC method, 
one calculates all accessible states equally populated in 



the decay of thermalized system into N fragments. The 
microscopic thermodynamics used here describes the de- 
pendence of the volume of 6N-dimensional phase space 
on globally conserved quantities (energy, mass, charge, 
...) and external constraints (like the spatial volume) 
to be defined by the first stage of the reaction. Within 
the microcanonical ensemble method, an explicit treat- 
ment of the fragment positions in the occupied spa- 
tial volume allows for a direct extension of the MMMC 
code [Q to the case of non-spherical shapes. Here, the 
source deformation is considered as an additional ex- 
ternal constraint. Main results of our paper will be 
given for the source described as an axial ellipsoid : 
{x/R x f + {y/R y f + {z/R z f = 1, with R x = R y £ R z . 
We assume that the freeze-out density of deformed sys- 
tem is the same as that of a spherical system with the 



radius R sys = (RxRyRz 
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i.e., the volume of deformed 



system is conserved. This condition changes neither the 
pass scheme nor the weight w r due to the accessible vol- 
ume of the fragments in the Metropolis scheme of calcula- 
tions (!]]. On the other hand, it means that the ellipsoidal 
source shape depends on one additional parameter : the 
ratio of ellipsoid axes 1Z = R x /R z . The ratio 1Z < 1 cor- 
responds to the prolate form, while for the oblate form 
one has 7Z > 1. 

An essential feature of non-spherical systems is that 
the deformation 'costs' some extra energy E^ef which is 
proportional to a change of nuclear surface with respect 
to the spherical shape. Since we do not consider shape 
evolution of the system but rather the influence of source 
shape on its thermodynamics, this energy Edef will be 
inaccessible for thermal motion and may be disregarded 
in the total energy balance. However this point should 
be kept in mind if one tries to refer to the real values of 
energy pumped into the system. 

The source deformation will noticeably affect the mo- 
ment of inertia and, together with the Coulomb energy 
which is calculated exactly for every multifragment con- 
figuration of non-spherical nucleus, becomes very impor- 
tant for describing rotating systems. As to the general 
scheme to account for the total angular momentum and 
the calculation of the statistical weight w p i of the config- 
uration in the rotating frame, we closely follow Ref. |q] 
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to be realized in the available code, [j For each spatial 
configuration of fragments, part of the total energy goes 
into rotation and hence the temperature of the system 
will slightly fluctuate. We take into account fluctuations 
of the moment of inertia arising from fluctuations in the 
positions of fragments and light particles. In the cal- 
culation of statistical decay of fragmenting system, the 
angular velocity of the source is added to the thermal 
velocity of each fragment. 

In calculating all accessible states within the standard 
MMMC method, the source should be averaged with re- 
spect to the spatial orientation of its axes, which is as- 
sumed to be homogeneously distributed in the whole 4-7T - 
solid angle jjj . If the angular momentum of fragmenting 
system, caused by the dynamical first step of the reac- 
tion, is strictly conserved then not all these states will 
be accessible and a formal averaging over the 4tt - solid 
angle will result in the violation of angular momentum 
conservation |0J|]. In the considered reactions, the an- 
gular momentum vector is perpendicular to the reaction 
plane. So, we disentangle in the MMMC code the beam 
direction (the z - axis) and the rotation axis (the x - axis). 
The rotation energy is then : L 2 /2J X = L x /2J X , where 
J x is the rigid-body moment of inertia with respect to the 
x axis. Averaging over the polar angle 9 is not consis- 
tent with the angular momentum conservation. On the 
contrary, averaging over 2-7T in the angle <f> corresponds 
to averaging over azimuthal angle of the reaction impact 
parameter and should be included. Averaging over ro- 
tation angle ip around L depends on the considered re- 
action, namely on the relationship between a rotation 
time : r ro t = J x /L Xl and a characteristic life-time of the 
source r c . For a system with high angular momentum 
when T rot <C t c , the full averaging in < if> < 1~k should 
be performed. In the opposite limit when r rot ^ t c , only 
states with ip ~ are accessible. Below we shall consider 
both these limiting cases. 

In the HI collisions, a part of the total energy can be 
stored in the compression energy of pre-formed source 
which during the collective (isentropic) expansion is 
transformed into the kinetic energy of fragments. In a 
strict thermodynamic sense, such an expanding system 
in not in equilibrium. However, in the quasi-static expan- 
sion, i.e., when the time scale involved in the expansion 
is larger compared to the equilibration time, the system 
may be considered to be infinitesimally close to the ther- 
mal equilibrium and consistently treated likewise an equi- 
librated system under the action of a negative external 
pressure whose magnitude is equal to the flow pressure. 
Such an approach has been applied for describing the flow 
effect in multifragmentation within a quantum-statistical 
model giving an estimate of about 20 MeV/A for a 



*In our version of the MMMC code the program error made 
in Q is corrected (see also Ref. Q]). 



maximal applicable flow energy. Consistent treatment of 
this effect within the microscopic approach would require 
an introduction of proper weight factors in the MMMC 
code. Such a work is now in progress flb]| . However, to 
get some insight into the influence of collective flow on 
the multifragmentation process, we shall mimic this effect 
by simply adding the blast velocity to the thermal ve- 
locity of each particle/fragment for any event simulated 
by the Metropolis method. This approximate procedure 
does not affect the Metropolis pass scheme in the origi- 
nal code. Assuming v ~ r, a simplified scaling solution 
of the non-relativistic hydrodynamic equations describ- 
ing the radial expansion of a spherical source provides 
the following general form for the radial velocity profile 

I3S 

v h {r) = v (r/R ) a . (1) 

Here vq and Rq are strength and scale parameters of the 
flow respectively, and the power-low profile is character- 
ized by the exponent a , which commonly is taken in the 
interval : 0.5 < a < 2. For the non-spherical expansion, 
the velocity profile may have a more complicated form 
and depends on the direction. But, even in the case of 
axially symmetric expansion, the scaling relation ([!]) was 
successfully applied for describing the profile of the trans- 
verse expansion velocity (lj|ll| . Below it is assumed that 
the radial expansion at the freeze-out point is described 
by eq. (|l|) . In the hydrodynamic interpretation, the scal- 
ing parameter Rq corresponds to the size of the system 
at the initial time of the scaling regime and, therefore, 
it should be less than the effective radius of the spheri- 
cal source at the freeze-out point, i.e., Rq < R sys . The 
strength parameter vq is then equal to the blast velocity 
at the surface of this effective sphere. As can be seen 
from eq. (|l|), the kinetic energy of fragments produced 
in the 'interior' will be mainly sensitive to the variation 
of the exponent a while the total collective flow energy 
may be fitted by different combination of all these three 
parameters. Note that the average collective energy of 
expansion, similarly as the deformation energy, is not in- 
cluded in the value of the total excitation energy. 

As an example, let us consider the multifragmenta- 
tion of 197 An having the angular momentum L = 40h 
and the excitation energy 6 A ■ MeV. These parame- 
ters correspond to the source formed in central Xe + Sn 
collisions at 50 A ■ MeV |l4|]. All calculations were car- 
ried out at the standard break-up density p po/6, 
what gives R sys — 12.8 fm for the radius of spherical 
197 Au nucleus. We consider two ellipsoidal shapes char- 
acterized by the axis ratio 7Z — 0.6 (prolate shape) and 
1Z = 1/0.6 = 1.667 (oblate shape). We have found that 
none of the observables related to the fragment-size dis- 
tribution is sensitive to the deformation of fragmenting 
source at such high excitation energies. The c.m.s. an- 
gular distribution of the largest fragment , i.e., the frag- 
ment with the largest charge (Z = Z max ) , is shown in 
Fig.l. In the absence of collective expansion, the angu- 
lar distribution is isotropic for oblate configurations and 
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has small forward - backward peaks for prolate config- 
urations if the averaging around L is performed in the 
whole available interval < -0 < 2ir. For the 'frozen' 
spatial configuration (ip — 0), the 'deformation effect' is 
clearly seen: for the prolate form there are strong forward 
- backward peaks , while for the oblate form the heaviest 
fragment is predominantly emitted in the sideward direc- 
tion (9 cm = 7r/2), like in the scenario of hydrodynamic 
splashing. This deformation effect is definitely not due 
to the angular momentum as can be seen by comparing 
Figs, la, lb (L = 0) with Figs, lc, Id (L = 40%). The 
collective expansion (a — 2) enhances the deformation 
effect. One may notice a strong enhancement of forward 
and backward peaks in the prolate case and the appear- 
ance of a strong peak at 8 cm = tt/2 in the oblate case. 
Similar features can be seen also in the cumulative angu- 
lar distributions of all IMF's but the relative amplitude 
of the deformation effect in that case is smaller. 
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FIG. 1. Angular distribution of Z max in the c.m.s. Plots 
on the l.h.s. correspond to averaging over the whole available 
interval < ip < 2-k, whereas plots on the r.h.s. correspond 
to the 'frozen' configuration ip = 0. (a),(b) : L — 0, v b — 0; 
(c),(d) : L = 40R, v b = 0; (e),(f) : L = 40h, v b = 0.08c, 
a=2 . 

Large sensitivity to the source shape is expected in the 
analysis using global variables on an event-by-event basis 
|l5j . Here we shall associate the global variables with the 
momentum tensor : 



(2) 



where p[ is the ith Cartesian coordinate [i = 1, 2, 3) of 
the c.m.s. momentum p( v > of the fragment v. The sum 
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FIG. 2. distribution. The plots on the l.h.s. cor- 

respond to the averaging over the whole available interval 
< < 27r, whereas the plots on the r.h.s. correspond to 
the 'frozen' configuration ip = 0. (a),(b) : L — 0, vu = 0; 
(c),(d) : L = 40ft, v b = 0; (e),(f) : L = 40h, v b = 0.08c, 
a=2. 

in (|) is running over all IMF's (Z > 3). The factor 
depends on the physical interpretation which one wants 
to give to the tensor (||). We use 7 = l/2m(„), where 
mi v \ is the mass of fragment v. The tensor Qij can be 
represented as an ellipsoid in the momentum space. The 
shape of this ellipsoid can be described by three axes 
and its orientation can be fixed by three angles in the 
3D-momentum space. This is usually done by referring 
to the eigenvalues < Ai < A2 < A3 (Ai + A 2 + A 3 = 1) 
of the tensor and to the Euler angles defining the 
eigenvectors ei,e2,e3. From various possible combina- 
tion of these parameters defining global variables |lf| 
, we consider here the sphericity : s = (3/2) (1 — A3), 
the coplanarity : c = (v / 3/2)(A2 — Ai), the aplanarity : 
a — (3/2) Ai, and the flow angle 0/z ot0 defined as an an- 
gle between ei and the z-direction (the beam direction) 
in the c.m.s. 

Distribution over the flow angle is presented in Fig. 
2. Different physical situations considered in Fig. 2, are 
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exactly the same as in Fig. 1. Even though &fi ow char- 
acterizes now all IMF's rather than the most sensitive 
largest fragment, the difference in the source shape man- 
ifests itself already for a vanishing blast velocity = 0. 
The whole effect is extremely sensitive to the presence of 
collective expansion (see Figs. 2e and 2f) and is enhanced 
furthermore for the 'frozen' configuration (ijj = 0). The 
maximal size of genuine angular momentum effects can 
be seen by comparing Figs. 2b [L = 0) and 2d (L = 407i). 
In general one expects that the cos(Ofi ow ) - distribution 
for highly central events is uniform 14 . However this 
naive expectation may be altered by many effects, such 
as the spatial shape of the source, the non- vanishing col- 
lective radial expansion , the high angular momentum 
and last but not least the detection bias. 





FIG. 3. (a) Sphericity , (b) coplanarity and (c) apla- 
narity distributions for L = 40ft. Bold lines correspond to 
Vfo = 0.08c, q=2. The continuous lines show prolate 7^=0.6 
shape whereas dashed lines show oblate 7?.=1.667 shape. 

Fig. 3 shows the distributions over sphericity dP/ds, 
coplanarity dP/dc and aplanarity dP/da for different 
source deformations. By construction, these distribu- 
tions do not depend on -0-averaging. The curves plotted 
with bold lines correspond to non- vanishing collective ex- 
pansion which reveals different ellipsoidal source shapes. 
The reason for this sensitivity can be seen from (Q) . De- 



pending on the source shape, a different number of frag- 
ments can be placed in the two regions : r < Rq and 
r > Ro (we put everywhere Rq = 0.7R sys ), in which the 
expansion acts differently. For the same reason, these dis- 
tributions are insensitive to the radial expansion for the 
spherical source. Particularly interesting are the spheric- 
ity and coplanarity distributions where the evolution of 
distributions with is clearly different for prolate and 
oblate source shapes. Lack of sensitivity of the aplanarity 
distribution to the deformation effects and to the expan- 
sion happens accidentally for chosen parameters. One 
should also mention that the sensitivity of dP/ds and 
dP/dc distributions to the source deformation and their 
insensitivity to the time-scales involved (the tp - averag- 
ing) provides interesting and supplementary informations 
to those contained in the Z maa; -angular distribution and 
in the cos(Ofi ow ) - distribution. 

The Z-dependence of average kinetic energy E% of 
IMF's for diffferent source shapes and different a- 
parameters is shown in Fig. 4. The calculations have 
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FIG. 4. Kinetic energy of fragments is plotted as a function 
of Z. (a) Prolate shape, L — 40ft. Different curves correspond 
to : Vfo — (circles), and : = 0.08c, a =1/2 (diamonds), 
1 (triangles), 2 (squares), (b) The same as in (a) but for the 
oblate shape, (c) The same as in (a) but for the spherical 
shape, (d) Comparison between the experimental data and 
the calculation (a=2) for both prolate and oblate freeze-out 
configurations and the blast velocities : = 0.088c (prolate 
shape), Vfo — 0.063c (oblate shape), chosen as to reproduce 
the experimental mean kinetic energy per nucleon of IMF's : 
5.2 ± O.lMeV/nucl. 

been done for the 'frozen' configuration (^=0). The col- 
lective expansion energy is a dominant contribution to 
the average kinetic energy of fragments, and its value 
differs noticeably for prolate and oblate source shapes. 
It is of interest to note that the average kinetic energy 
of fragments exhibits a flattening or even a maximum at 
large Z (see , e.g., the curve for a = 2). The precise 
position of this maximum depends also on the deforma- 
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tion (see Figs. 4a and 4b). An attempt to reduce this 
large observed kinetic energy to angular momentum ef- 
fects results in a value of the average angular momentum 
(L w 640S.) which is completely unrealistic for selected 
central events 14 . Fig. 4c shows Ek(Z) for the spheri- 
cal source shape. Fig. 4d compares results of the present 
model with the experimental data for central Xe+Sn col- 
lisions at 50 A-MeV jlj]. All calculated events have been 
filtered with the INDRA software replica, and then se- 
lected with the experimental centrality condition : com- 
plete events (i.e., more than 80% of the total charge and 
momentum is detected) and Qfi ow > 7r/3. 

In conclusion, external constraints on the shape of 
cquilibrized fragmenting source have been considered 
within the extended MMMC method. Due to the change 
in the Coulomb energy for deformed freeze-out configu- 
ration, the shape effect is clearly seen in the IMF's an- 
gular distributions (Z max - angular distribution) as well 
as in the fi ow - distribution. A surprising interplay be- 
tween effects of non-spherical freeze-out shapes and the 
memory effects of nonequilibrium phase of the reaction, 
such as the rotation and the collective expansion of the 
source, has been demonstrated for the first time. The 
influence of shape on rotational properties of the system 
is not only reduced to the modification of the momenta 
of inertia. The limits on the averaging interval over the 
angle -0 about the rotation axis, which are defined by 
the appropriate time scales, affect strongly the angular 
observables and are able to enhance strongly the 'shape 
effect'. These constraints may be important for certain 
observables used in experimental procedures of selecting 
specific class of events. Other striking finding is that 
the collective expansion allows to disclosure the source 
shape in the analysis using global variables as well as in 
the study of Z-dependence of the average kinetic energy. 
The latter observable is independent of ip- averaging but, 
unfortunately, it is specified by a poorly known profile 
function ([!]). Nevertheless, the careful analysis with eq. 
(|l|) might shed some light on the problem how different 
fragments are situated in the freeze-out configuration. 

In the experimental analysis, if the average kinetic en- 
ergy of fragments is fixed by an appropriate choice of 
vjj, Rq for each source deformation , then the shape of 
Ek(Z) contains information about the exponent a in the 
parametrization (1) and to the lesser extend about the 
deformation of the source. The form of the dP/ds and 
dPj dc distributions permits then to find the average de- 
formation of the fragmenting source. Having fixed the 
parametrization (1) and the deformation of the source, 
the analysis of the angular distribution of Z max and/or 
the fi ow distribution gives an access to the informa- 
tion about the time-constraints (the limits on the tp- 
averaging) in the multifragmentation process. We have 
discussed here the freeze-out shape effect for only one 
value of the excitation energy. Certainly, the manifesta- 
tion of this effect in observables is energy/ angular mo- 
mentum dependent and the modifications of the MMMC 
method, presented in this Letter, open a promising way 



to follow up the shape evolution of equilibrized source 
in the broad range of incident energies. Alongside with 
the observables discussed here, it would be interesting to 
study the velocity correlations between fragments which 
are sensitive to the source shape at the freeze-out. Such 
a work is now in progress Eof . 
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